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ABSTRACT 

We count the supersymmetric bound states of many distinct BPS monopoles 
in N = 4 Yang-Mills theories and in pure N = 2 Yang-Mills theories. The 
novelty here is that we work in generic Coulombic vacua where more than 
one adjoint Higgs fields are turned on. The number of purely magnetic 
bound states is again found to be consistent with the electromagnetic du- 
ality of the N = 4 SU(n) theory, as expected. We also count dyons of 
generic electric charges, which correspond to 1/4 BPS dyons in = 4 
theories and 1/2 BPS dyons in = 2 theories. Surprisingly, the degen- 
eracy of dyons is typically much larger than would be accounted for by a 
single supermultiplet of appropriate angular momentum, implying many 
supermutiplets of the same charge and the same mass. 
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1 Introduction 



The purpose of this paper is to count BPS monopoles and dyons in supersymmetric 
Yang-Mills theories. Such computations have been performed by many authors using 
moduli space dynamics of monopoles |TJ, ^ |3|, f|, |5|] . Recently |jL |7|, 0, however, it was 
realized that the monopole dynamics in a generic vacuum is qualitatively different 
from the old moduli space dynamics of Manton |§ employed in most such endeavors, 
where the low energy dynamics of monopoles were considered when only one adjoint 
Higgs field is turned on, while supersymmetric Yang-Mills field theories come with 
2 or 6 such scalars. This restriction disallows static interaction between monopoles 



10], so that all interaction comes from nontrivial coefficients of kinetic terms. In a 
generic vacuum with more than one adjoint Higgs turned on, monopoles of the same 
type still have no static force among them, but dynamics of monopoles of distinct 
type could have a static potential. In this article, we solve various index problems as 
a first step towards counting all BPS states. 

In the old monopole dynamics, the Lagrangian one finds is a pure sigma model 
with the moduli space as the target manifold. Classically, one solves for geodesic 
trajectories to find classical orbits of monopoles. Quantum mechanically, the Hamil- 
tonian is a square of a supercharge which can be regarded as a Dirac operator acting 
either on the spinor bundle or on the Clifford bundle over the moduli space, 

Q „ -z 7 - Vm . (1) 

Supersymmetric bound states, for example, would be found as normalizable spinors 
or forms on the moduli space that are also zero modes of this Dirac operator. 

The new supersymmetric low energy dynamics is obtained by augmenting old 
moduli space dynamics with a set of supersymmetric potential terms, and was written 
explicitly in Ref. ]7|, §. In principle, the question of BPS states must be reconsidered 
in the new dynamics. (One example of states that cannot be probed in the old 
formalism is the now well-known 1/4 BPS states [11, [T^, |I3||.) In this new setting, the 



supercharges of the low energy dynamics will again be interpreted as Dirac operators 
on the moduli space, which is now twisted by some triholomorphic Killing vector 
field, say K, 

Q~l m {-iV m -K m ). (2) 



The corresponding bosonic potential is precisely half the squared norm of K pA 



We count the number of normalizable states annihilated by a Dirac operator, weighted 
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by ±1 for the chiral and the antichiral states respectively. The resulting integer is the 
index of the Dirac operator. In fact this index would be infinite or ill defined without 
further restriction of the domain of the Dirac operator. We restrict the problem to 
each charge eigensector before counting the zero modes of the Dirac operator. Thus, 
effectively, we will be computing an equivariant index with L 2 condition. This gives 
information on the existence and the degeneracy of dyonic bound states for each 
electric charge sector. 

Understanding monopole dynamics in generic vacua is particularly significant in 
the context of N = 2 Yang-Mills theories, because many (1/2) BPS dyons exist only 
when both adjoint Higgs are turned on. By ignoring the potential term in the low 
energy dynamics, one would in effect be searching for a bound state in the vacua 
where the state cannot exist as a supersymmetric one-particle state. In the language 



of Seiberg-Witten fL5| , one would be looking for it on one side of a marginal stability 
domain wall in the vacuum moduli space, while the bound state in question exists 
only on the other side of the domain wall. For example, most dyons that become 
massless at special hypersurfaces in the Seiberg-Witten moduli space, are of this type 



Tq , [T7|, |T8| and cannot be probed by old moduli space dynamics. 

In Section 2, we review the supersymmetric quantum mechanics with potential 
for the case of 4 and 8 supercharges. We isolate various involutions, with respect 
to which the index is defined. Section 3 introduces the explicit form of the moduli 
space metric that governs the dynamics of BPS monopoles. The only known moduli 
space for arbitrarily many monopoles is the case of all distinct monopoles, and this 
is the case for which we will compute the index explicitly. Section 4 recalls a recent 
explicit computation of two-monopole bound states in SU(3) theories and presents 
the resulting value of the indices. In section 5, we finally delve into the computation 
of the index by using a Fredholm deformation of the Dirac operator in question. The 
computation reduces to that of a superharmonic oscillator in 4 dimensions, whose 
index is computed explicitly. Section 6 translates the results to degeneracies of various 
dyonic and purely magnetic bound states and checks its consistency with anticipated 
nonperturbative physics. We close with a summary. 



2 Supersymmetric Sigma Model with Potential 

In this section, we briefly review the supersymmetric sigma-model quantum mechan- 
ics with potentials. These quantum mechanics have been introduced as the low energy 
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dynamics of monopoles in pure N = 2 Yang-Mills field theory and in N = 4 Yang- 
Mills field theory |7], [8]. They also appeared in other systems such as the dynamics 
of instanton solitons ||19| . 



2.1 Quantum Mechanics with 4 Real SUSY 

The SUSY dynamics we consider is a sigma model with potential, whose Lagrangian 
is written compactly as 

C = X - [g mn z m z n + ig mn X m D t X n - g mn G m G n - zV m G n A m A n ) , (3) 

where D t X m = A m + r^,i n A p . The target manifold must be hyperKahler, which 
means that there are three covariantly constant complex structures J( s > satisfying 
quaternionic algebra, 

j(s)j(t) = _ S st + e stuj(u)^ ( 4 ) 

and the Killing vector field G should be triholomorphic. 

C G g = 0, C G J {S) = 0. (5) 

Introducing vielbein and defining X E = X m e E % which commute with all bosonic 
variables, the canonical commutators are 



{X E ,X F } = 5 EF . (6) 

We can realize this algebra on spinors on the moduli space by letting X E = j E /V2, 
where •y E are gamma matrices. (Since the moduli space is hyperKahler an equivalent 
quantization is obtained using holomorphic differential forms.) The supercovariant 
momentum operator, defined by 

7T m = p m — -^U mEF [X , X F ], (7) 

where u} 1 F E is the spin connection, then becomes the covariant derivative acting on 
spinors 7r m = — zV m . Note that 

[vr m , A n ] = «r™ p X p , 

["'mi 7I~n] ~^R"mnpqX X . (8) 
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The supersymmetry charges take the form 

Q = A m (7r m — G m ), 
QW = A m J«> n -G n ), (9) 

which obey 

{Q,Q} = 2(H-Z), 
{Q {s \Q {t) } = 26 st (H-Z), 
{Q,Q (s) } = 0. (10) 

Here the Hamiltonian TC and the central charge Z are given by 

H = \(^Ti m Vgg mn v n + G m G m + iV m G n \ m \ n y (11) 
Z = G m -n m - % -(V m G n )\ m \ n . (12) 



Note that the operator iZ is the Lie derivative Cq acting on spinors (see e.g., P0| ) 



C G = G m V m + iv m G n [ 7 "\ 71- (13) 

o 

The SUSY quantum mechanics comes with a natural Z 2 grading defined by the op- 
erator, 

r 2 = n2 1/2 A s = Il7 E , (14) 
which anticommutes with the Dirac operator, 

D = V2Q = 7 m HV m - G m ). (15) 

This pair defines the Witten index that counts the difference between the number of 
bosonic states and the number of fermionic states annihilated by the supercharge. In 
fact, the index is defined in each superselection sector with fixed Z, and effectively 
counts the difference in the numbers of BPS states of given central charges. The 
index will be denoted collectively by X 2 . See Section 5 for detailed computation of 

x 2 . 

2.2 Quantum Mechanics with 4 Complex SUSY 

When the number of supercharges and the number of fermions double, we obtain the 
following form of sigma model with potential, 

-g m nGTG r l - iV rn G In ^(Q^r), (16) 
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where %p m is a two component Majorana spinor, 7 = 02, 7 1 = «0"i,7 2 = — ias, 



tp = ip T 'j . The operator fij's are defined respectively by f2 4 = 6™ 7^3, f2s = <5™7, 2 



and f2 s = iJ^ m n5ai3 for s = 1,2,3. The supersymmetry algebra again requires 
the manifold to be hyperKahler. As in the previous subsection, the G 7 's must be 
triholomorphic Killing vector fields. 

When quantized, the spinors ip E = e^tp"" 1 with vielbein e,^, commute with all the 
bosonic dynamical variables, especially with p's that are canonical momenta of the 
coordinate z's. The remaining fundamental commutation relations are 



[z m , Pn ]=i5™, 

{^,^} = S EF 5 a ,. (17) 
Define supercovariant momenta by 

TTm = Pm ~ ~U EF m ljj E J°7p F , (18) 

where wef™, is the spin connection. The N=4 SUSY generators in real spinors can 
be written as, 

Qa = OW. - (7°^^) m ^, (19) 

g w = (j(*V)> m - (^j^n^rcl . (20) 

These charges satisfy the N = 4 complex superalgebra: 

Qp} = {Q { :\ Qp ] } = 25 aP U - 2( 7 °7 1 ) a/3 Z A - 2( 7 °7 2 W Z s , (21) 

7J f — * lap Zs, IVq'i^ f — "la? 

)f} = 2^Z 1 , {Qf^} = 2 1 l, 
where 7i is the Hamiltonian, and the Zi's are central charges 



{Q a , Qf } = 2 7 ^ Z s , {QS\Qf>} = 2 7 ^ Z 3 , (22) 

{gi 2) ,g? ) } = 2 7 ° zx, {g<, 3 \g£ ) } = 2 7 ° J z 2 , (23) 



= G^iXm - l -V m Gi ^ m 7 V n - (24) 

The sigma- model without the potential possesses an 5*0(5) R-symmetry which is 
explicitly broken by the G lJ s. The G / 's transform as 5 of SO(5)r. 

The complex form of the supercharges is often useful. To this end, we introduce 
ip m = ^(V'r — ^T) an d define Q = -j^iQi — iQ2)- The supercharges in ( fL9|) can be 
rewritten as 



Q = tfr* n - ^ m {Gt - iGl) - i G s m (J {s) ^ m , (25) 

s=l 

gt = <p™* m - ^{Gl + iG 5 m ) + i £ G s m {J {s) vT • (26) 
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The charges and are analogously defined from (|20|). The positive definite 
nature of the Hamiltonian can be seen easily in the anticommutator 

{Q,Q j } = {Q {s \Q^} = 2H, (27) 

while the central charges appear in other parts of the superalgebra. For instance, we 
have 

{Q,Q} = -Z A + iZ 5 , 
{Q\Q ] } = -Z i -iZs. (28) 

Once we adopt this complex notation, it is natural to introduce an equivalent ge- 
ometrical notation. Defining the vacuum state |0) to be annihilated by yj* m 's, and 
using the 1-1 correspondence, 

(^ m V m2 --^ mfc )|0) «-> dz mi A dz m ' 2 A • • • A dz mk , (29) 

we can reinterpret (p m as the exterior product with dz m , and ip* m = g m n { P* n as the 
contraction with d/dz m . The supercharge operators can be rewritten as, 

Q = id LQ4_iQ5 + ^ /Ql\ ~\~ ^j(2)(£.2) "I" ^J(3)((j3)> 

Qt _ j^t _ t G4+ . G5 - Uj(l) (G l) - !ij(2)( G 2) - iij(3)( G 3), (30) 

where is the contraction with the vector field if, and its conjugate t^- is the exterior 
product by the 1-form obtained from K by lowering its indices. 

The SUSY quantum mechanics admit a canonical Z 2 grading, which in the geo- 
metrical notation of (129|) is defined on fc-forms by 

r 4 ee (-l) k . (31) 

or equivalently by 

The involution T4 anticommutes with all supercharges and determines the usual Wit- 
ten index, X4. 

In some special limits, however, there could be an additional Z 2 grading. Suppose 
that we have only one nonzero G 1 , say G 5 . The operators 



T ± = H(Vi<p E ±V=i<p* E ). (33) 
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anticommutes with the Dirac operators defined as, 

At = %Q ± Q ] = {%^ m ± ^ m )(vr m t G 5 J = d - l g s ± i(S - 4), (34) 
the square of which is 

D 2 ± = ±2i(H =f Z 5 ). (35) 

So the Z 2 gradings define an analog of the signature index for each choice of sign 
and for each charge-eigensector. A given state with nonzero Z 5 can be annihilated 
by one of D± at most, and in fact must break at least half of the supercharges. The 
corresponding indices will be denoted by if. 

In Section 5, we will compute both X4 and Xf in such a special limit with only one 
of five G /, s present, which we can take to be G 5 without loss of generality. For I4, 
we may take any one of D± as the Dirac operator, since T4 anticommutes with both. 
A standard index theorem will then allow us to deduce I4 in more general setting. 



3 Moduli Spaces 

Moduli space dynamics of monopoles decompose into the interacting relative part and 
the non-interacting "center of mass" part. The latter corresponds to a 4-dimensional 
flat metric of the form, 

g cm = AdX 2 + Bde T , (36) 

where A is a three-vector. Since we are interested in establishing existence of bound 
states, this part of the dynamics will be ignored for the most part. 

The free center-of-mass sector generates two kinds of quantum numbers, never- 
theless. One is the overall, conserved U(l) charge, and the other is a supermultiplet 
structure generated by the fermionic partners of A and £r- The resulting degen- 
eracies, 4 and 16 for N = 4 real and complex supersymmetric quantum mechanics 
respectively, correspond to the smallest possible BPS multiplet of the underlying 
SUSY Yang-Mills field theories with 8 and 16 supercharges, respectively. 



3.1 Distinct Monopoles 

A simple case of this dynamics involves a collection of distinct monopoles in SU(n) 
gauge theories. The interacting part of the moduli space metric is a simple generaliza- 
tion of four-dimensional Taub-NUT metric |2l| . Without loss of generality, consider 
a collection of k + 1 distinct monopoles, whose magnetic charges are given by an 
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irreducible (sub)set of simple roots, (3i, . . . , Pk+i- The simple roots satisfy relations 
(31 = 1, p a ■ j3 a+ i = —1/2, and f3 a ■ j3 a+ b = for b > 1. The relative part of the 
corresponding metric is 

4tt 2 

g = C ab dr a - dr b + —(C )ab{di) a + cos 9 a d4> a ){dipb + cos9 b d(p b ), (37) 

where the matrix C for the relative moduli space is[] 

2tt 1 

C a b = ^ab H n^ab — • (38) 

e 2 r a 

The 3-vector f a is the relative position between the a th and (a + l) th monopoles, 

r a = %a+l %ai (39) 

while the angles ip a of period Arc are related to the U(l) phases of each monopole, 
£ a 's (of period 2tt), by 

^ d d d 

dlp a d^a+1 d£ a 

/fe+1 \ Q k+1 Q 

(y m \° = Ym a ^-. (40) 



where £t is a coordinate that appears in free center-of-mass part of the dynamics and 
m a is the mass of the a th monopole. 

For a generic reduced mass matrix fi, the triholomorphic Killing vector fields of 
this geometry are exhausted by 

Ka = (41) 

so the vector fields G and G 1 are linear combinations of K a 's with constant coeffi- 
cients; 

G = e y a c K c , 

c 

G 1 = eya{K c . (42) 

c 

The electric charges are measured by the charge operators, 

~iC Ka , (43) 



lr The coupling constant e will be assumed to be positive without loss of generality. 
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whose (half-) integer eigenvalues will be denoted by q a . In terms of the simple roots 
[3 a , the electric charge of a dyonic state with charge g a 's is 

e(+qi + q 2 + q3 + --- + qk + ra/2)0i+ 

e(-gi + q2 + <?3 H \~ qk + n/2)/3 2 + 

e(-qi-q2 + q3 + --- + qk + n/2)(3 3 + 

e(-qi -q 2 -q3 + --- + q k + n/2)(3 4 + I 44 ) 

e(-qx - g 2 - ?3 qk + n/2)(3 k+1 . 

where the integer n comes from quantization of an overall U(l) angle and should be 
even or odd when 2 J2a qa is even or odd, respectively. 

3.2 Unit Noncommutative Instanton 

A simple deformation of the above moduli space appeared in another context recently, 
where one considers low energy dynamics of an instanton soliton in the 5-dimensional 
noncommutative U(k + 1) Yang-Mills theory |22]. This happens because an instanton 



in S 1 x R 3 can be regarded as a collection of k+1 distinct monopoles of the underlying 
Yang-Mills theory |23|. When we compactify the theory on a circle of radius R, the 



nontrivial part of the moduli space of a single instanton soliton is given by the metric 



El 



g = g2 (C ab dr a ■ df b + (C'^abidtpa + cos 6 a d<j) a ) (di/) b + cosd b d<f) b )^ , (45) 

where e is the 5-dimensional Yang-Mills coupling. The matrix C for the relative 
moduli space is 

1 1 

Cab = v ab + 5 ab h =r— , (46) 



r„ 



Efa-2irC/R\ 



where ( encodes the noncommutativity. The matrix v is determined by the Wilson 
line along 5* 1 that breaks the gauge symmetry to U(l) n . When the supersymme- 
try of the underlying field theory is maximal with 16 supercharges, the low energy 
dynamics of the instanton is given by our SUSY quantum mechanics with 4 com- 
plex supercharges. When the field theory comes with 8 supercharges, the instanton 
dynamics is described by the SUSY quantum mechanics with 4 real supercharges. 



4 Bound States of Two Distinct Monopoles 

For a pair of two distinct and interacting monopoles, the dynamics have been solved 
for supersymmetric ground states in each charge eigensector. The geometry is that 
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of a Taub-NUT manifold which comes with a single triholomorphic Killing vector 
field K\. Accordingly, there is only one conserved U(l) charge, qi, which labels 
superselection sectors. 

In the pure N = 2 Yang-Mills case, define a\ = 47r 2 ai/e 3 /x where // is the reduced 
mass and ai is defined by G = ea\K\. The normalizable wavefunctions had been 
constructed by Pope in another context [23|, and the number of dyonic bound states 
of charge q\ was found to be [Q 

2|<7i|, (47) 

if < q± < a or a\ < q\ < 0, and 

0, (48) 

otherwise. For each q±, the solutions belong to the same chirality spinors, and thus 
contribute to the Witten index equally. Thus the Witten index X 2 in each charge 
eigensectors are 



To 




< |<7i| < | Cfci | and < qdi 
otherwise 



(49) 



For a pair of distinct monopoles in N = 4 Yang-Mills f25| , the five G 7 's must be 
proportional to the single triholomorphic vector field Ki = d/dipi- We may rotate 
them into a single triholomorphic vector field, say G I=5 = ea\K\, upon which we can 
define a± similarly as above, = 4n 2 ai/e 3 fi. The degeneracy is found to be 



1, 



for purely magnetic state (gi = 0), while for dyons 



(50) 



(51) 



when < \q±\ < \ai\, and zero otherwise. All solutions are self-dual differential 
forms, when we take the convention that the curvature tensor of the moduli space is 
self- dual. 

When the central charge Z 5 = ea\q\ of the state is positive (negative), the bound 
state is annihilated by D + (£L) only, while for Z$ = 0, the state is annihilated by 
both. For given qi, we find 



11. 



1 gi = 

|gi| < \q\\ < |5x| and qia\ > 

otherwise 



(52) 
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and 



t: = l 



1 qi = 

\qi\ < \q±\ < \di\ and qid\ < 

otherwise 



> . 



(53) 



The Witten index X4 counts the number of even forms minus the number of the odd 
forms. Of solutions with q\ 7^ 0, half are even and the other half are odd, so we find 
that 

'1 qi = 

J 4 = } (54) 

,0 Ql + 

regardless of a±. 



5 Index Computation 

We would like to put a lower bound on the number of bound states in the above 
SUSY dynamics by computing indices. The index problems can be quite involved, 
given that the quantum mechanics involve many degrees of freedom with complicated 
interaction terms. However, the problem can be simplified by utilizing the invariance 
of the index under certain deformations. In this section we will use the invariance of 
the index under Fredholm deformation to simplify our index computations. Before 
proceeding with the computation, however, we need to restrict to the regime where 
a massgap exists. 



5.1 Massgap 

When restricted to specific charge eigensectors, the operators above may exhibit 
two drastically different behavior. For small charges, the sector has a massgap; the 
continuum part of the spectrum is bounded below by a positive gap. For large 
charges, the massgap disappears. This is the reason why there is an upper bound on 
the electric charge q\ of bound states of two monopoles. In the two-body problems, 
the condition for the massgap to exist in a sector with electric charge q 1 is 



0i < 



47r 2 I ail 



(55) 



where the bosonic potential is generated by a single triholomorphic vector field G = 
ea\K\. When we consider many distinct monopoles, the condition for the massgap 
to exist is equally simple: 

|<7c| < \a c \, (56) 
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where 

471-2 k 

a c = r )cba b , (57) 

e 6=1 

with G = ea c K c . In the quantum mechanics with four complex supersymmetries, the 
same holds true provided that only one G 1 , say, G 5 = ea c K c is turned on. We will 
compute the indices, X 2 , X 4 , if assuming that all of these conditions hold.f] 

5.2 Index Generalities 

First we recall basic definitions. 

Definition 5.1 A bounded linear operator L : E\ — > E 2 between two Hilbert spaces 
is Fredholm if there exists a bounded operator P : E 2 —>■ E\ such that PL — J 1; and 
LP — I 2 are compact operators. Here Ij denotes the identity map on Ej. 

The operator P in the above definition is called a parametrix. We will be interested 
in the case where L is a Dirac operator. In this case, although Dirac operators are 
unbounded on L 2 , we may trivially make L bounded by taking E\ to be the closure 
of the domain of L with respect to the norm (graph norm) 

H/lliU^ imi 2 + !|£/!l 2 , 

where unsubscripted norms denote L 2 norms. 

If L is a Dirac operator on a compact manifold, then it is well known to be 
Fredholm. In the compact case, one takes, for example, P to be the Green's operator, 
Q defined to be the unique operator satisfying: 

(i) Q annihilates the kernel of U . 

(ii) The range of Q is orthogonal to the kernel of L. 

(iii) LQf = f for / e the image of L. 

Then Q is bounded by 1 + A^ 1 ^ 2 , where Ax is the first nonzero eigenvalue of L^L. 
Also, QL — I\ — 111, and LQ = I 2 — II2, where TEi and II2 denote the orthogonal 
projections onto the kernels of L and 1/ respectively (and are finite rank and thus 
compact operators). Hence P = Q satisfies all the conditions of the definition. 
2 If five G 1 = e ^ c a^KcS are involved, the massgap condition generalizes to 

( 9c ) 2 <]T(^) 2 , (58) 

7=1 

where arc defined similarly as above for each G 1 . However, the Indices Tf- are not well-defined 
unless all G^'s are proportional to each other. We will discuss such generic cases in Section 6. 
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In the case of a Dirac operator on a noncompact manifold the preceding construc- 
tion of a Greens operator may fail to yield Fredholmness for several reasons. The 
kernel of L or L* may fail to be finite dimensional, making one of the projections 
not a compact operator. Also, if there is no gap in the spectrum, Q will fail to be 
bounded. These deficits are all avoided, however, under the assumption that the 
essential spectrum of L^L is bounded away from zero. (We recall that the essential 
spectrum includes the continuous spectrum and any eigenvalue of infinite multiplic- 
ity.) Then the kernels are finite dimensional and Q is again bounded by 1 + X 1 1 ^ 2 , 
where Ai is the smallest nonzero element of the spectrum of L^L. It is well known 
that the essential spectrum is bounded away from zero whenever L^L has the form 
A + V, where for two positive constants c and R, V satisfies V(x) > c for x outside 
a fixed compact set. All the operators we consider in this paper have this form. 



A basic result in index theory (eg |26| p. 122), is the following. 



Proposition 5.2 Let L t , t e [0, 1] be a continuous family of Fredholm operators. 
Then index(L Q ) = index(Li) . 

Thus one can sometimes deform an index problem to a more tractable index 
computation. To avoid potential confusion, we recall the notion of continuity assumed 
in the above proposition. L t is a continuous family of operators if for each s and for 
every e > 0, there exists 5 > so that \\L t F — L S F\\/\\F\\ < e for all nonzero F 
if \t — s\ < 8. In particular, we note that we require the 5 to be independent of 
F. Hence, for example, the super harmonic oscillator in one variable ip\4- -\- fax, 
cannot be continuously deformed to ipi 4- + tp 2 by scaling away the interaction term. 
If one uses the graph norm for ipi-£; + ip2, then the oscillator is unbounded and hence 
clearly cannot be deformed to a bounded operator. If one instead uses the graph 
norm for the oscillator, it is easy to see that all frequencies give equivalent norms and 
by construction, multiplication by x (as a map to E 2 is continuous in each of these 
norms. Hence the deformation 

L t =ip 1 ^- + ^ 2 ({l-t)x + t) 

is continuous. The limit operator, however, is not Fredholm as a map from E\ — ► E 2 
even though it is easy to show that it is Fredholm if the oscillator graph norm is 
replaced by the ipi + ^2 graph norm. 

In analyzing continuous families of operators L t it is often useful to utilize also 
families of parametrices P t . If, however, we choose P t to be the Greens operator Q t 
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of L t then we will be plagued by the possibility that if eigenvalues converge to zero, 
Qt will become unbounded. Hence, for no other reason than to avoid such problems 
of bounding P t , it is useful to define a modified Greens operator 

g Lt , £ = g t (i-n Lt , £ ), 

where Hi t ,e denotes the projection onto the A < e eigenspaces of L[L t . This operator 
is bounded by 1 + e -1 / 2 and is a parametrix as long as e lies below the essential 
spectrum. 

5.3 Deforming the Index 

For several of the index problems we will be considering, it seems likely that one 
can simply deform the given operator into a standard superharmonic oscillator and 
then immediately deduce the index. There are some minor issues fitting such a 
deformation into a continuous family. We will not treat those here because for one of 
our index computations - that of the noncommutative instanton - there is no single 
model operator to which to deform. Instead we will use the deformation invariance 
of the index to localize all the problems to an elementary computation around the 
zeros of our triholomorphic vector field G. 

The case of interest to us then is D a Dirac operator anticommuting with an 
involution r, L the restriction of D to the +1 eigenspace of r, and E x and E 2 the spaces 
of sections of the associated bundles with finite graph and L 2 norms respectively. 
In this context, Fredholmness follows from the conditions in the preceding sections 
guaranteeing a mass gap (i.e., bounding the essential spectrum of D^D away from 0.) 

The deformations we will consider involve replacing G by TG for some T large. 
This fits into the above framework without modification since T > 1 ensures preser- 
vation of the mass gap. Moreover, scaling G is clearly continuous because the norm 
of G is a bounded in the given metric. Recall that G enters the Dirac operators in 
the form of operators, 

\ m G m (59) 

or 

(v^f iv^DG™ (60) 

which are Clifford multiplications by G. Denote these operators by G. We see that 
the sup norm of the difference between two Dirac operators (associated to TG and 
SG) is bounded by \\{TG - SG)f\\ < |(T - S)\ x \G\ sup x ||/||, which clearly gives 
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the desired inequalities for the continuity of the deformation. We note that even had 
the metric allowed for unbounded \G\, we still would have G bounded as an operator 
from Ei (equipped with the graph norm) to E 2 , as in the oscillator example of the 
previous section. 

In addition, we will modify the metric on certain compact subsets. This modifi- 
cation may change the actual domain and range of our operator. For example r, and 
hence its eigenspaces may vary with the metric. Nonetheless, we may choose quasi- 
isometries between them. Thus if we have Fredholm operators Dt ■ E\(T) — > E 2 (T) 
and quasiisometries hi(T) : Ei — > Ei(T) then the index of h 2 {T)~ 1 DTh\(T) is T 
independent by the proposition and is equal to the index of Dt since the index is 
unchanged under composition with bounded operators with bounded inverse. We 
note, although we will not need it here, that the condition that hi be quasiisometries 
may be relaxed to the condition that the eigenvalues of hi and grow at most 
polynomially (subexponentially even) in distance from some choice of origin. This 
is an easy consequence of the fact that the Fredholm estimate implies exponential 



decay of the elements in the L 2 kernel of Dt- (See p7|). As we will use these decay 
properties, let us recall them in a crude form now. 

Suppose we have N points y^i = 1, • • • , N and a Hamiltonian of the form H = 
A + 4T 2 V, with V(x) > 1 if \x - y t \ > 1, i = 1, • • • , N. Suppose also that Hf = X 2 f 
for some small constant Ao, and / G L 2 , say with L 2 norm 1. Let 

\x\ m := mmi<i<jv|a; - 

Then e (T2 ~ A o )1/2|l ' |m / e L 2 , and the L 2 norm of | e (T2 ~ A o) 1/2 l x 'l m /| restricted to the 
exterior of the balls of radius R > 1 about the y\ is finite and bounded by 4e^ T2 ~ x ^ 1/2R . 



(This is not sharp. See [57J for sharper statements.) In particular, we observe that 
the L 2 norm of / restricted to the complement of the balls of radius 2R about the yi 
satisfies 

||/|^|| 2 <4e-^ 2 - A o) 1/2 « (61) 

Hence, / is concentrated near the zeroes of V. 

Let Dt denote our Dirac operator with G replaced by 2TG and the metric modified 
to be Euclidean in a ball of radius 10R some R » 1 about each zero of \G\ 2 . This 
metric modification allows us to compare Dt to a model Dirac operator which agrees 
with Dt near the zeros of G and has known index. Assume, as we may by replacing 
G initially by a suitable multiple, that D t Dt has the form A + 4T 2 V, with V{x) > 1 
in the complement of the balls of radius 1 about each zero of G. As in the previous 
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section, for e below the continuous spectrum of D T , Hd Tj£ denotes the projection onto 
the A < e eigenspaces of D T . Then 

Index/A^ = J dx tr rHT,e{x,x). (62) 

Using (|6l| ) we see that for d e := rankno Ti£ , 

/ dx tr tU Dt)E (x, x) = f dx tr rn DT . e (a;, x) + 0(4e" 2c(T2 " e)1/2 ), (63) 

J J\x\ m <2c 

for some choice of c. Hence it suffices to bound d e independently of T large and to 
compute the integral of tr rH^x, a;) over \x\ m < 2c in the large T limit. 

First we estimate d e . Let D T f = A 2 /, for A < e, and ||/||l 2 = 1. Let Qt denote 
the Green's operator for the super harmonic oscillator (SHO) which agrees with D T 
in a neighborhood of radius 4R about the zeros of G. Let p R denote a cutoff function 
supported on a ball of radius 2R where D T is the SHO and identically one on a ball 
of radius R. Then (p R f — Qt-DtPrJ) is in the kernel of the SHO. Denote its norm 
by a and introduce a unit vector v in the kernel of SHO such that 

(pnf ~ QTDlpRf) = av. 

Observe that Qt has sup norm < T _1 . Now consider the equality 

D 2 T p R f = X 2 PR f+[A,p R ]f. 

By our assumptions, the right hand side is 0(A 2 ) + 0(e~ TR ) (not sharp). Hence 
\\Pnf ~ av \\ — 0(X 2 ) for A > 0(e~ TR ^ 2 ). Setting, for example, e = 1/T, we have then 

\\f-av\\=0(l/T 2 ). (64) 

Moreover, such an inequality is true for any vector in the image of Hd t .i/t- AVe 
conclude then that rank Hd t ,i/t is no larger than the dimension of the kernel of 
the SHO (times the number of zeros of G). This bounds dyx and completes our 
demonstration that it suffices to compute the trace over a bounded region. 

In the following, for simplicity of notation we will consider the case of a single 
zero for G, but the general case follows similarly with only notational complications. 

Let St and te denote the Euclidean Dirac operator and involution which agree 
with Dj* and r near the zeros of G. Let Ft denote the Greens operators for S^- To 
define a parametrix for Dt, introduce Qt, the Greens operator for Dt, and let Pt be 
the modified Greens operator 

P T ■= Q T (J - U Dt ,i/t)- 
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Define 

h := Index( J D+) - Index(S£) = Tr([D+P T - S^F T ] - [P T D+ - F T S$]). 

Then we wish to show that the integer 1\ = 0. By (|63|) , the above traces can be 
approximated for R, T large as 

h = Trp R ([D+P T - S+F T ] - [P T D+ - F T S+]) + 0(d 1/T e- RT/2 ). 

On the support of p R , = S£, hence we have 

h = Trp R [D+, P T - F T ] + 0{d 1/T e- RT ' 2 ). 

Using the cyclic property of the trace, we rewrite the first term on the right hand 
side of the above formula as 

Trp R [D+,P T - F T ] = -Tr[D + , p R ](P T - F T ) + J dx ViV{x)\ 

where V(x) is the vector with 

Vi(x) := Tr 7i (p i? ,(x)(P r - F T )(x,x)). 

The integral vanishes by Stoke's theorem, leaving 

h = - J ti[D + ,p R ]{x){P T -F T ){x,x)dx. 

We estimate the last term by converting it back into an expression involving the 
exponentially decaying projection operators. Write 

[D+,p R ](P T -F T ) 
= [D^p R ](F T S^ + II 5 + l/T ){P T - F T ) 

= [D+,p r ]Ft(S+Pt - S+F T ) + [D+,p R ](U s + )1/T )(P T - F T ). (65) 

The last term is 0(Te~ RT ) because 

[D^p R ](U s ^ 1/T ) = 0(e~ RT ) 

by ( |B~ID and because (Pt — Ft) has sup norm < T by construction. We separate the 
first term into two additional terms 

[D T i Pr]Ft(S t Pt — StF t ) 
= [D+, p R }F T {D+P T - S+F T ) + [D + , p R }F T (S} - D+)P T . (66) 
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The first term is again exponentially decreasing because (D^Pt — S^Ft) is a difference 
of exponentially decaying projection operators and Ft is uniformly bounded. We can 
compute Ft explicitly and it is 0(e~ TS ), where 5 is the distance between the support 
of [Dt,Pr] and the support of (S^ — Dt)- Hence all the terms are exponentially 
decreasing, and we deduce that I\ is exponentially decreasing. On the other hand, it 
is the difference between two integers and must therefore vanish. 
We summarize our results. 

Index(£#) = Index(^). 

When there is more than one zero of G, a minor variation of the same argument 
yields 

Index(£#) = ^Index(S' T (i) + ), 

i 

where St (2) is the local model for Dt at the i th zero of G. In order to extend our 
argument to this case, we must replace expressions of the form FtPt and Tl s + i/ t Pt 
in the previous expression by FtPutPt and II 5 + 1 / T PtitPt for some large n, because 
Ft and U s + X j T need not extend naturally to the full moduli space in the many zero 
case. This will introduce new error terms of the form V PnTPrPTtobeestimated. Our 
decay estimates can once again be used to show these terms are also exponentially 
decaying. 

5.4 Computing the Deformed Index 

We now use the deformation arguments of the preceding section to complete the index 
computations. We consider first the case of the quantum mechanics with 4 complex 
supersymmetries on a moduli space of dimension 4k and compute Xf. The other cases 
are very similar and follow with minor modifications. In the case, we have reduced 
the problem to computing the index of the operator S\/ e := d— Lq+i{S — t G ) acting on 
selfdual forms (r + = 1) on C 2k . Separating variables, we see that the index of Si is the 
product of the indices of the D c , c — 1, • • • , k, where where D c := d— i ac K c +i{d^ — 4 c x c ) 
(no sum over c) acting on selfdual forms on C 2 . Using the deformation invariance of 
the index again, we may assume a c = 2. 

This latter index is easy to calculate exactly as follows. We compute 

D C D\ + D\D C 

= A + \2K C \ 2 - {d, t\ Kc } + i{t 2Kc: d} - i{d\ i\ Kc } - {d\ t 2Kc }. (67) 
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Let z\ and z 2 denote complex coordinates on C 2 . Then 2K C = iJ2 2 =i(z c -^. — z c~§t\ 
Hence, \2K C \ 2 = |^i| 2 + \z 2 \ 2 . In the coordinate frame, on (p,q) forms we have 

i{i2K c , d} = -i{d\ t\ Kc } = iC 2Kc = -(p -q)~ 2K c /i, 

and 

{d, i 2K } = idzidzi + idz 2 dz 2 , 

{d\ i>2K c } = idz\dz\ + idz 2 dz 2 . (68) 

Hence, we have 

D C D\ + D\D C 

= A + \z\ 2 - 2{p - g) - 4^ C A - {d, l\ Kc } - {dt, i2Xc }. (69) 
The functions 

/(a, 6, c, d) := (9 21 - - Zl )\d Z2 - z 2 ) c {d- Z2 - z 2 ) d e-^ 2 / 2 , (70) 

with a, b,c,d> 0, span the eigenspace of A+|z| 2 with the eigenvalue 2(a+b+c+d)+A. 
We compute commutators to obtain 

2K c f(a, b, c, d) = i(—a + b — c + d)f(a, b, c, d). 

Therefore, on the algebraic span of f(a, b, c, d) we have 

D C D\ + D\D C = 4(a + c) + 4 - 2(p - g) - {d, i\ Kc } - {d\ i 2Kc }. (71) 

This vanishes if and only if a = c = and the differential form coefficient of f(a, b, c, d) 
takes one of the following forms: 

dzi A dz 2 , (72) 

or 

1 + [idz\ A dz,\ + id^2 A dz 2 )/2 — dz 1 A d^i A dz 2 A dz 2 , 

or 

(dzi + idz 2 A dz 2 A dzi/2), 

or 

(d^2 + idz\ A d^i A dz 2 /2). 
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We have, therefore, an infinite dimensional kernel to our operator before taking 
into account the constraint on charge. We now recall that we wish to restrict to the 
space 

2q c = C 2K Ji = {p-q)+ 2K c /i. 

With the above normalization of K c , g c 's are integers or half-integers. On /(0, b, 0, d) 
this imposes the constraint 

2g c = (p-q) + (b + d). 

The index of D c is thus given by the number of ways to choose nonnegative integers 
b and d so that 

b + d = 2q c - (p - q) 

with p — q = 2 or plus twice the number of ways to choose nonnegative integers b 
and d so that 

b + d = 2q c -l. 

There are 8q c such solutions for positive q c , one such solution for q c = 0, and none 
for negative q c . All of these solutions are self-dual, so the index of D c is 

8q c q c > 

q c = > . 
Qc < , 

Note that this result assumes a positive coefficient of K c . For a negative coefficient, 
the computation proceeds exactly as above, provided that we make the following 
exchanges of coordinates, 



Index (D c ) 



1 





z 2 <-> z 2 



(73) 



This maps K c to — K c , and flips the sign of q c in the charge constraint above. In 
other words, the sign condition in the index formula is really on a c q c for each j. Thus 
the index is 

a c q c > 



Index (D c ) 



8\q c 
1 




a c q c = 
a c q c < 



for each c. 

Thus, whenever there exist a massgap, the index Xf is 

a c q c > 
a c q c = 
a c q c < t 
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where the sum is over the zeros of the potential. Note that the index is nonvanishing 
only if all a c q c (no summation) are nonnegative. The states in the kernel of the Dirac 
operator must be annihilated by Ti — Z as well, and the central charge Z of the states 

e]Ta c g c >0 (74) 

c 

equals the energy. 

Computation of X~, appropriate for those states with positive central charge, 
proceeds similarly. In fact, this problem can be mapped to that of J+ by 



<p 


* 

-> up 




-> V 




- -K c 



(75) 

The net effect is to flip the sign condition on the charges q c , so 

(f 8|g c | a c q c < 

n 1 « c ^=o 
c { a c q c > t 

whenever a massgap exists. The sum is over zeros of the potential. The energy of 
the contributing states is — e^ c a c q c > 0. 

We consider next the same set of operators but now restricted to the +1 eigenspace 
of T4. We see that the terms with p + q even are in the +1 eigenspace of T4, and the 
terms with p + q odd are in the —1 eigenspace of T4. This leads to zero index for all 
nonzero q c . When q c = 0, we get a solution with p — q = 0. There is only one of 
these. The index of r 4 is then 

:;!}) 

where the sum is over the zeros of the potential. 

Finally, we consider the minor modifications necessary to compute I 2 . Once again 
a separation of variables allows us to reduce the index of the Euclidean operator to 
a product of indices of operators B c , I — 1, • • • , k on C 2 . In coordinates, B c has the 
form 

2 d d 

acting on the +1 eigenspace of 4A1A2A3A4. This choice of the Dirac operator corre- 
sponds to positive coefficients a c — 2 for all c = 1, . . . , k. Then in a covariant constant 
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frame, 

2B C B\ + 2B\B C = A + \z\ 2 + AiK - A 2j _iA 2i . 

i 

Acting on the algebraic span of f(a, b, c, <i), 

2B c £t + 2B\B C = 4(a + c) + 4 + 4* £ A^_iA 2j -. 

i 

This has infinite dimensional kernel spanned by the product of /(0, 6, 0, d) and a co- 
variant constant spinor in the (one dimensional) intersection of the —1/2 eigenspaces 
of iAiA 2 and iA 3 A 4 . 

The charge constraint in a covariant constant frame takes the form 

2q c = 2K c /i — iAiA 2 — 1X3X4. 

Acting on the above basis elements of the kernel of 2B C B\ + 2B\B C this reduces to 

2q c = b + d+l. 

Counting as before this yields a 2q c dimensional kernel which lies entirely in the +1 
eigenspace of 4AiA 2 A 3 A 4 . Hence, the index of B c is 



2\q c \ a c q c > 

a c q c < 



Index (B c ) - 

Thus, whenever there exist a massgap, the index X 2 is 

^=e(ii 



2\q c \ a c q c > 

a c q c < 



Finally, we conclude the index computations by noting that they reproduce the 
four- dimensional results summarized in the previous section. In fact, the wavefunc- 



tions found in Ref. [25| and in Ref. 24| can be seen easily to reduce to the superhar- 



monic oscillator wavefunctions above in the limit of a 3> q. 

6 BPS Bound States 

The above index computations count differences in the number of ground states with 
respect to Z 2 involutions r±, T4,t 2 , 

Index = n + — ri- (76) 
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where n± are the number of ground states with r eigenvalue ±1. We are actually 
interested in the sum + instead, for which one needs a more refined understand- 
ing of the dynamics. For the case of r 2 and r±, we anticipate n_ vanishes by itself. 
Such a vanishing theorem is shown rigorously for the simplest cases in Appendix. We 
will assume in this section that n_ = holds true for r 2 and t± in all cases, and 
compare the results to what are expected on physical grounds. 

6.1 N = 4 Yang-Mills Theories 

The supersymmetric quantum mechanics with four complex supercharges describe 
dynamics of monopoles in N = 4 Yang-Mills theories. Recent studies of D-branes 
indicates the following three possibilities for dyonic bound states of monopoles. 

• The state is 1/2 BPS in the Yang-Mills field theory. These states would be 
annihilated by all supercharges of the low energy monopole dynamics, which is 
possible only if the central charges in the relative part of the dynamics is absent. 
This is guaranteed when all relative electric charge g a 's vanish. In particular, 
this includes purely magnetic bound states. 

• The state is 1/4 BPS in the Yang-Mills field theory. These states would be 
annihilated by half of the supercharges of the low energy monopole dynamics 
and not by the other half. This is possible only if at least one central charge is 
nonzero. 

• The state is non-BPS. 

The index computation of the previous section tells us something about 1/2 BPS and 
1/4 BPS states, where we counted indices Z 4 and Xf in the special limit where only 
one G 1 , say G 5 , is turned on. Equivalently, we considered vacua where two Higgs 
fields are turned on. 

Of the three indices, only J 4 is robust against turning on more than one G T, s. The 
Dirac operator iQ ± would no longer anticommute with r± but does anticommutes 
with r 4 . Only Z 4 is a well-defined index in such generic vacua. Turning on additional 
G 1 always increases the massgap, and is a Fredholm deformation that preserves X 4 . 
Thus our result shows that, in generic vacua, 

X 4 = 1, (77) 
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when q a = 0, and zero otherwise. Since the central charge of the state that contributes 
to the index is zero, the state must be annihilated by all supercharges of the quantum 
mechanics and is a 1/2 BPS in N = 4 Yang-Mills theory.^ This is consistent with 
the existence of a unique magnetic 1/2 BPS bound state of monopoles in generic 
Coulomb vacua, which is expected from the SL(2, Z) electromagnetic duality. One of 
the generators of SL(2, Z) maps massive charged vector multiplets to purely magnetic 
bound states in 1-1 fashion. After taking into account the automatic degeneracy 16 
from the free center-of-mass fermions, the total degeneracy of these bound states 
is alway 16, which fits the N = 4 vector multiplet nicely. This purely magnetic 
bound state was previously constructed by Gibbons || in special vacua where all 
G^s vanishes. 

Existence of 1/4 BPS states are more sensitive to the vacuum choice and the 



electric charges. The existence criteria were first found by Bergman where 
he constructed these dyons as string webs ending on D3-branes. The first necessary 
condition is that the string web should be planar, which is equivalent to the condition 
that, effectively, only one linearly independent G 1 is present. This allows us to assume 
without loss of generality that only G 5 is turned on, as far as counting 1/4 BPS 
states are concerned. Thus, the computation of Xf in the previous section is directly 
applicable. 

Secondly, at each junction of the string web, the string tensions must balance 
against each other, which in the present language of low energy dynamics translates 
to the condition that the effective potential in the charge-eigensector is nonrepulsive 
along all asymptotic directions fT2]| ; 

|<?d < \a c \. (78) 

This second condition may indicate the existence of a minimal energy bound state, 
however, does not guarantee that the state would preserve some supersymmetry. 

Finally, a minimal energy configuration is supersymmetric when the orientation 
of string segments are consistent with each other. Say, if one fundamental string 
segment is directed to one particular direction, then another fundamental string in 
the same web must be directed the same way. The second string can point toward 



3 One might think that existence of this bound state is obvious since the potentials are all attrac- 
tive and also there exists a classical BPS monopole of the same magnetic charge. However, none of 
these guarantee the existence of BPS bound state at quantum level. In fact, the same set of facts 
are true for a pair of distinct monopoles in TV = 2 SU (3) Yang-Mills theory but we know that such 
a purely magnetic bound state does not exist as a BPS state B. 
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the opposite direction and still balance the tension, but such a combination breaks 
all supersymmetry. This orientation condition on the string web, is nothing but the 
condition that the product a c g c 's (no summation) are all of same sign. See figure 1. 
Thus, a 1/4 BPS dyon may exist only when \q c \ < \a c \ for all c and a c q c are all of 
same sign, at least one of which is nonzero. 




(0,1) 



(0,1) 



(0,1) 



(1,0) 

> — 



(-1,0) 

— < — 



(-1,0) 

— < — 



(a) (b) (c) 

Figure 1: Simple webs of (q,p) strings that represent dyons in N = 4 SU(4) theory. 
The filled circles represent D3-branes, while vertical lines are D-strings, and horizon- 
tal lines are fundamental strings. Configurations (a) and (c) preserve a quarter of 
supersymmetry that was left unbroken by D3's, while (b) breaks all supersymmetry. 
When translated to Yang-Mills field theory on D3's, the horizontal separations be- 
tween D3-branes are encoded in a c 's while the electromagnetic charges are determined 
by which string ends on which D3. 

The indices If- were computed with the massgap condition \q c \ < \a c \ to begin 
with, and yielded nonzero value only when all a c q c were of the same sign; For positive 
signs of a c q c , X+ ^ 0, while for negative a c g c 's, we have T~ ^ 0. The result is 
clearly consistent with the existence criteria set by the string-web construction, and 
furthermore gives us extra information beyond the string web picture. The index 
indicates that the degeneracy of such a 1/4 BPS state is 

16 x []Max{8|g c |,l}. (79) 

c 

The factor 16 arises from the free center-of-mass fermions. 

In the two monopole bound states, the number 8|g| is accounted for by four angular 
momentum multiplets of j = \q\, \q\ — 1/2, \q\ — 1/2, \q\ — 1 The top angular 



4 The first three suffices for \q\ = 1/2. 
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momentum \q\ in the relative part of the wavefunction has a well-known classical 
origin: when an electrically charged particle moves around a magnetic object, the 
conserved angular momentum is shifted by a factor of eg /Air. While fermions can 
and do contribute, the number of fermions scales with the number of monopoles, and 
not with the charge q a . In fact, it is most likely that the top angular momentum of 
such a dyonic bound state wavefunction is 

jtop = J2\lc\, (80) 

c 

for large charges, so that the highest spin of the dyon would be 

i+j top = i + £kd, (si) 

c 

after taking into account the universal vector multiplet structure from the free center- 
of-mass part. On the other hand, a 1/4 BPS supermultiplet with the highest spin 
jtop + 1 has the total degeneracy of 

16x8]T|g c |, (82) 

c 

which is much less than the number of states we found above unless all but one q c 
vanishes. Thus, this implies that there are many 1/4 BPS, thus degenerate, super- 
multiplets of dyons for a given set of electromagnetic charges. This is probably the 
least understood of our results. While one would expect to find degenerate states 
within a supermultiplet, there is no natural symmetry that accounts for the existence 
of many supermultiplets of the same electromagnetic charges and of the same energy. 

For large electric charges q a , the number of dyon supermultiplets scales as, at 
least, 

(j[Max{8\q c \A}j / (s^lq^ . (83) 

Proliferation of dyonic states of a given charge was anticipated by Kol some time 

ago in the context of string webs in 5 dimensions ||28|| . Because our computation 

was performed for a collection of distinct monopoles, which put some constraint on 

the possible topology of the related string web, it is not immediately clear to us 

whether we can make any sensible statement in the regime where Kol's prediction 

is applicable.^ Nevertheless, it is tantalizing that we found the number of states 

Kol anticipated exponentially large numbers of states, in fact, which is much more than our 
power like result. 
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increasing much faster than would have been expected from supersymmetry alone. 
It is not clear to us why this happens and what interpretation this may have in the 
Yang-Mills field theory. 

In the regime where \q c \ > \a c \ for some q c , we cannot rely on the current index 
computation. On the other hand, since even a single repulsive direction, i.e., \q c \ > \a c \ 
for some c, prohibits a bound state (supersymmetric or not), the unresolved question 
boils down to the marginal case, where \q c \ equals \a c \ for some c's while the others 
satisfy \q c \ < \a c \. The only state that must exist for sure is the purely magnetic 
bound state (q a = 0), which was constructed by Gibbons when a c = and which is 
necessary for SL(2, Z) invariance. The explicit construction of two-monopole bound 
states in Ref. seem to indicate that no dyonic bound state may form along such 
marginal directions, but this remains to be seen for multi-monopole cases. 

6.2 N = 2 Yang-Mills Theories 

In N = 2 Yang-Mills theories, a state could be either BPS or non-BPS. There is no 
such thing as a 1/4 BPS state. Dyons that would have been 1/4 BPS when embedded 
in N = 4 theories, are realized as either 1/2 BPS or non-BPS depending on the sign 
of the electric charges. The index of this Dirac operator was nonzero only when 

< q c < a c for all c 
or > q c > a c for all c (84) 

which gives us a possible criterion for BPS dyon to exist. This condition is similar to 
the condition for BPS dyons or monopoles to exist in N = 4 Yang-Mills theories but 
differs in two aspects. The first is that given a set of a c 's, all of which are positive 
(negative), the electric charge g c 's must be all positive (negative). 

The second difference from N = 4 case is that a purely magnetic bound state of 
monopoles does not seem to exist as a BPS state, even though there exists a classical 
BPS solution of such a charge. This feature was noted previously in Ref. |J, where 
bound states of a pair of distinct monopoles were counted explicitly. In fact, the 
index indicates that all relative q a must be nonvanishing for a BPS state to exist. 
Assuming the vanishing theorem, the number of BPS dyonic bound state under the 
above condition is 

4xl]2|g c |, (85) 

c 

The overall factor 4 is from the quantization of the free center-of-mass fermions. 
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For large electric charges we again observe the proliferation of supermultiplets. 
The top angular momentum and thus the size of the largest supermultiplet can grow 
only linearly with YJ \Qc\ which means that the number of supermultiplets with the 
same electric charges scales at least as 

(n 2 i?ci)/( 2 £i?ci) (86) 

for large g c 's. Again it is not clear to us what implication this has in the Yang-Mills 
field theory. 

In the regime where \q c \ > \a c \ for some q c , again we cannot rely on the current 
index computation. For the same reason as in N = 4 Yang-Mills theory, no bound 
state can exist if even a single repulsive direction (\q c \ > \a c \ for some c) exists, 
so the unresolved question boils down again to the marginal case, where \q c \ equals 
\a c \ for some c's while the others satisfy \q c \ < \a c \. Extrapolating from the explicit 
construction of two-monopole bound states in Ref. ||, we suspect that no bound 
state may form along such marginal directions. 

6.3 Ground States of a Noncommutative Instanton Soliton 

Supersymmetric ground states and excited BPS states of an instanton soliton in 
S 1 x R 3+1 can be counted similarly as above. The only difference as far as comput- 
ing the index goes, is that the potential has many zeros. For a single instanton in 
noncommutative U(n) theory, there are precisely n zeros of the bosonic potential, 
and near each of these zeros, the Dirac operator can be deformed to that of a su- 
perharmonic oscillator. One crucial difference in interpreting the result in physical 
terms comes from identification of the conserved charges. Of n possible conserved 
£7(1) charges, n — 1 relative charges are again electric charges. However, the overall 
conserved U(l) does not correspond to a gauge symmetry, and comes from translation 
of the instanton along S 1 . This last U(l) charge is just the Kaluza-Klein momentum 
along S 1 . 

Of particular interest are the quantum ground state of the instanton with no U(l) 
charges excited. In the maximally supersymmetric U(n) Yang-Mills theory, the index 
tells us that there are n distinct BPS supermultiplets of ground states. This result 
was anticipated in Ref. [jnj. With half as much supersymmetry, however, the index 
is consistent with no supersymmetric quantum ground state exist at all. 

Quantum states of instanton soliton in R 4 was previously studied in the com- 



mutative setting |29], |30|]. In particular, the absence of a quantum ground state of 
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instanton soliton in the nonmaximal supersymmetric Yang-Mills theories, has been 
observed from string-web construction. 

7 Conclusion 

By computing indices and assuming vanishing theorems, we counted supersymmetric 
bound states of arbitrarily many distinct monopoles in N = 2 pure Yang-Mills the- 
ories and also in N = 4 Yang-Mills theories. The relevant low energy dynamics are 
supersymmetric sigma-models with potential(s), where the supercharges preserved 
by supersymmetric bound states can be interpreted as Dirac operators twisted by 
triholomorphic Killing vector fields. An obvious generalization of this computation 
is to include hypermultiplets in N = 2 Yang-Mills theories, but it goes beyond the 
scope of this paper. 

Counting of 1/2 BPS states in N = 4 Yang-Mills yielded a result consistent with 
electromagnetic duality of the theory. In particular, the necessary purely magnetic 
bound states of distinct monopoles are all accounted for in SU(n) theories. While 
this result is not surprising, it is still significant in that this was shown for the first 
time in all generic Coulomb vacua of the Yang-Mills theory. In contrast, distinct 
N = 2 monopoles do not seem to bind at all unless all possible relative charges are 
turned on. 

Existence criteria for 1/4 BPS states, previously found in the context of string- 
webs, are also faithfully reflected in the index formulae. On the other hand, the 
degeneracy of most 1/4 BPS dyons is shown to be much larger that one would have 
expected from a single 1/4 BPS supermultiplet with a physically reasonable angular 
momentum. N = 2 Dyons of the same electromagnetic charges as 1/4 BPS dyons 
of N = 4 theories, could be BPS or non-BPS, depending on the signs of the electric 
charge. We also counted the degeneracy of such N = 2 BPS dyons, which shows 
similar proliferation of supermultiplets. This phenomenon is not understood at the 
moment. It should be also interesting to see how the degeneracy behaves when both 
electric and magnetic charges are large. 

We are grateful to Jerome Gauntlett for stimulating discussions. We thank Aspen 
Center for Physics and also the organizers of the workshop, "The Geometry and 
Physics of Monopoles" , where this work was initiated. The work of M.S. is supported 
by NSF grant DMS-9870161. 
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Appendix 



Here, we prove a vanishing theorem for r 2 and r± on four-dimensional moduli space. 
Let us consider r 2 first. Because of the triholomorphic Killing conditions on G, dG 
is self-dual and does not couple to antichiral spinors. Then the Dirac operator is a 
simple Laplacian; 

- D m D m (87) 

when acting on antichiral spinors. Using the standard trick of sandwiching this 
operator by a hypothetical zero mode ^ and its complex conjugate we find 

= " [l[dz m V*V m D m V = - [l[dz m g mn (D n y)*(D m y) (88) 

m m 

where the possible boundary term vanishes by itself since the massgap forces \I/ to be 
exponentially small at large distances. Therefore, 

= D m V = (V m - iG m )^ (89) 

everywhere. This modified connection is still unitary as G m is real. Hence, ^ is 
covariant constant with respect to metric compatible connection and is therefore of 
constant norm. Such an / cannot be normalizable on an infinite volume space unless 
it is identically zero, which proves the vanishing theorem in four dimensions for r 2 . 

The case of t± can be handled similarly. Let us recall that differential forms can be 
thought of as a tensor product of two spinors. With an appropriate sign convention, 
we can identify various sectors of the former with those of the latter as follows 

• self dual even form — > [chiral spinor]® [chiral spinor] 

• self dual odd form — > [chiral spinor]® [antichiral spinor] 

• antiselfdual odd form — > [antichiral spinor]® [chiral spinor] 

• antiselfdual even form — > [antichiral spinor]® [antichiral spinor] 

On antiselfdual even form, then, the Clifford action of dG is trivial. Since the self- 
dual curvature does not couple to antiselfdual forms either, its action is also trivial. 
Thus, the square of the Dirac operator becomes a simple Laplacian again, 

D% = -D m D m (90) 

with D rn = V m ^iG 5 m . By the same logic as in the spinor case, therefore, no antiself- 
dual even- form solution can exist. Finally this also shows that antiselfdual odd- form 
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solution does not exist; Unless the central charge vanishes, a solution generates other 
solutions via the action of broken supercharges. The broken supercharges are the 
linear combination of Q and orthogonal to D±, so that it flips while preserving 
t±. Thus, the number of odd-form solutions equals the number of even-form solu- 
tions, in each r± eigensectors, whenever the central charge is nonzero. This proves 
the vanishing theorem for r± for sectors with nonzero U(l) charges q c . 
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